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Abstract 

Starting from full quantum field theory, various mean field approaches are 
derived systematically. With a full consideration of external source depen- 
dence, the stationary phase approximation of an action gives a nuclear mean 
field theory which includes quantum correlation effects (such as particle-hole 
or ladder diagram) in a simpler way than the Brueckner-Hartree-Fock ap- 
proach. Implementing further approximation, the result can be reduced to 
Hartree-Fock or Hartree approximation. The roll of the source dependence in 
a mean field theory is examined. 
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One of important problems in nuclear physics is understanding systematically the nuclear 
systems ranging from the structure of stable finite nuclei to a very hot and dense system 
which may occur in a high energy heavy ion collision or in a neutron star. A nuclear system 
is a strongly interacting many body system and thus the description of the system becomes 
complicated. A modern approach to the study of a nuclear system is based on the relativistic 
field theory in terms of relativistic nucleons interacting each other by exchanging mesons. 
The simplest approach of such a theory is the so-called relativistic mean field approximation 
(RMF) or quantum hadrodynamics (QHD) which describes a nuclear system in terms of 
nucleonic Dirac field interacting with classical meson fields |^. In this approach, a nuclear 
system is composed of relativistic nucleons whose self energy is determined through meson 
fields which are generated by the nuclear density. 

In spite of the success of this simple model in describing various properties of nuclear 
system with effective interaction [|]-|^, it is failed to reproduce nuclear saturation prop- 
erties from the free nucleon-nucleon interaction. The calculations yield correct binding 
energy with too large density or vice versa forming the Coester band |lTO|,p!I|. Even though 



the Dirac-Brueckner-Hartree-Fock (DBHF) approximation including ladder diagram [|T^-[T^ 



can achieve the nuclear saturation for a nuclear matter |T^JT^, it failed in reproducing the 



ground state properties for finite nuclei using the nucleon-nucleon interaction forming a 



Coester band [IS, IS]- The DBHF is essentially the RHF with correlation effects (ladder 
diagrams) determined through the Bethe-Salpeter (BS) equation for a nucleon-nucleon in- 
teraction. Two coupled self-consistent equations, Dyson's equation and BS equation, make 
the DBHF calculations very complicated for a finite nucleus |T^. Thus we need to search for 



a description that is simpler than the DBHF and can consider various quantum correlation 
effects. On the other hand, various calculations [^,20| show that the exchange term, the 



vacuum fluctuation, and the ladder diagrams can not be treated as perturbation with re- 
spect to the Hartree term or each other exhibiting non-perturbative characteristics of strong 
nuclear interaction. 

For a systematic study of various mean field approaches, it is necessary to investigate the 
basic assumptions underlying to each of the mean field theories starting from the quantum 
field theory. In this paper, we will consider full external source dependence of a stationary 
phase approximation (SPA) of the path integral for a nuclear system with arbitrary external 
sources. The added source terms bring the quantum correlations into a mean field approach 
in a simple way. Depending on the interaction considered, this approach can give a correction 
to the DBHF |^| and also reveals the origin of non-perturbative relationships among the 



various mean field approximations. 

Up to few hundreds MeV energy, a nuclear system may be described as a system of 
nucleons interacting each other via meson exchange. For a nuclear saturation property, we 
at least need to include a scalar meson for a long range attraction and a vector meson for a 
short range repulsion. The simplest Lagrangian describing a nuclear system can be written 
in the form of 

C{x) = ip{x) [ijf.d^' - M]iIj{x) + i[a^(^(a;)a^(/?(x) - m^^\x)] + gilj{x)Tip{x)^{x). (1) 

Here i/j is the nucleon field and Lp represents meson fields (scalar field and vector field here). 
The r is the unit matrix for a scalar meson and it is the Dirac matrix for a vector meson 
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field. More detail differences between scalar and vector meson fields are irrelevant for the 
discussions in the context of this paper. 

In quantum field theory, any physical observable can be represented as a function of field 
operators. Introducing arbitrarily small external sources as 



C{X,J) = C{x) — '0(x)J(x) — J{x)lp{x) — Jrrlx)ip{x), 

the expectation value of an operator Oitlj^^^Lp) can be obtained by 



W{J) \dJ' dJ' dJ, 



(2) 



(3) 



where the transition amplitude W{J') can be obtained through a path integral p^^4 

/ H{J)dtT^ 



W{J) 



(4) 



Here T represents time ordering and J stands for J, J, and/or J^. Notice here that J{x) is 
not the conjugate function of J{x)] both are independent arbitrary functions. If we use the 
single particle state representation, then the matrix element {C{x, J)) = \C{x, J') \ 
becomes a functional of the single particle wave functions and if{x). Now the quantum 
field theory for a nuclear system has been reduced to a problem of finding the corresponding 
transition amplitude W{J') or action S{J). 

For a Lagrangian of the form of Eqs.(|l|) and (0), we can integrate Eq.(^) over the meson 
field Lp. However the result becomes fourth order in the nucleon field if) and thus further 
integral is not feasible. On the other hand, we can integrate over the nucleon field first, 
but the resulting nonlinear equation of (p prevents further integration. Thus we need to use 
some approximate method. However, since we can not use a perturbative method due to 
the strong interaction, we are forced to use a classical or semi-classical trajectory as the 
simplest nonperturbative method pl|-p5|l. If we treat all the fields as classical ones, then 
JT")) in Eq.(^) becomes the Lagrangian £ of Eq.(0) with c-functions ipi^x), and 
J{x). To treat a nuclear ground state as a Slater determinant of occupied single nucleon 
levels, we should treat the nucleon field as a quantum field. Then the external sources J[x) 
and J(x) should also be treated as q-functions having the same quantum characteristics as 
the nucleon field with infinitesimal amplitudes to keep the source terms in the transition 
amplitude ||21[]. Similarly, J mix) should be a q-function if we consider a quantized meson 
field ^{x). 

The stationary phase approximation (SPA) or steepest descent method of the path inte- 
gral Eq.(^) or equivalently the variation of the Lagrangian with respect to each field gives 
the corresponding Euler-Lagrange equation for a classical trajectory as 



\i^^,d^ - M \ gV^{x)\^{x) = J{x), 
d^d^ + m^ if{x) = gi>{x)T^jj{x) - Jm{x). 



(5) 
(6) 



The equation of motion for is given as the conjugate of Eq.(|^). Using the meson propagator 
D{x — x'), we can eliminate meson fields from Eqs.(|^) and (^). The result becomes 
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t-ff,d^ -M + g'{ d^x'D{x - x')^{x')T^{x') T 



ip{x) 



= J{x) + gTi){x) j d^x'D{x - x')JUx'), (7) 

for a ground state of a nuclear system. This is a nonlinear equation of the nucleon field. 
Since the equations of i/j (Eq.(0)) and -0 are coupled nonlinear equations, their solutions 
would be functions of infinite order in the external source J', i.e., J, J and Jm- Usually this 
source dependence was neglected, and thus the mean field theory gave the RMF of Walecka 
|T|] or the Hartree-Fock (RHF) approximation only. Since Eq.(|^ should be satisfied order 
by order in J', we may expand the solution as a power series of JT"; 

ij{x) = u{x) + i)i{x) + ij2{x) + 4)3{x) + ■■■■ (8) 

Here U{x) = if^oix) is the source independent solution, which is non-zero only for occupied 
states, and ipn{x) is the n-th order term of the external sources J, J and Jm, i-e., i^nix) = 

The source independent solution U{x) can be obtained from Eq.(|^) by neglecting all the 
external source dependent terms; 



(i7^a^ - M) + / (^j c/VD(x - x')U{x')VU{x')^ V 



U{x) 



0. 



(9) 



Notice that the source independent field U{x) is the conjugate field of U{x). If we treat 
meson field Lp in Eq.(||) as a classical field, i.e., replace the right hand side of Eq.(|^) by the 
corresponding expectation value as it has been done by Walecka, then U {x')TU {x') inside 
the integral of Eq.(H) becomes a classical quantity (u{x')TU{x')'j and is independent of 
U{x) appearing outside the integral. This is the relativistic Hartree (RH) or the RMF of 
Walecka If we treat the meson field as a quantum field, then both the U{x') and U{x) 
are the same field and thus we have a nonlinear equation for U. To find the solution U{x) 
using an iterative method, we can linearize this equation as 



d'^x' A{x,x)U{x') 



0, 



{i-f^d^ -M)+g' d^x"D{x - x")U{x")TU{x")T 



S{x — x') 



- g^TU{x)D{x-x')U{x')T. 



(10) 



fill 



Here the minus sign of the last term in Eq.([Tl|) originates from the anti-commuting property 
of the nucleon field. Notice here that A{x, x') is hermitian and the self-consistent field 
U{x) in a nuclear system is different from the Dirac field in a free space. The U{x) and 
U{x") appearing in Eq.([l0|) through A{x,x') of Eq. ([TT|) are now treated as independent 
quantities from U{x') of Eq.(^) in the iterative calculation for the solution U{x). This 
is the relativistic Hartree-Fock (RHF) approximation. The second term of A{x, x') is the 
direct Hartree contribution to the nuclear self energy and the third term is the exchange 
Fock contribution to the self energy. In Ref. 0, it was shown that we can obtain RHF using 
the quantum operator algebra without using the concept of linearlization explicitly. 

On the other hand, from the n-th order terms of J in Eq.(0), we obtain coupled equations 
for Tpn and t/^^. Their solutions are 
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tpnix) = j d^x'B-\x,x')Jn{x') 

- J d^x'B-\x,x')TU{x') / d^X2 / d'^X3D{x' - Xs)Jn{x2)A-\x2,X3)TU{x3), (12) 

4'n{x) = J d'^x' Jn{x')B^^{x' , x) 

- J d^Xij d'^X3U{x2)TA~^{x2,X3)Jn{x3) / d'^x' D{x2 - x')U{x')VB-^{x' ,x), (13) 

where 

( fc """ _L I ^ J. I ^ J. n 

Jn{x) = -g'^ XI H H rV'ni(a;) / d'^x'Dix - x')lljn2{x')TllJn3{x') 
ni=0 n2=0 n-i=0 

+ J{x)5n,i + gTiJjn-i{x) J d'^x' D{x — x')Jmix'), (14) 

„ 71—1 71—1 71—1 

Jn{x) = -g"^ / d'^x'D{x' -X) ^ ^ Yl ^nA^')^'^n2{x')^ni{x)T 



n— 1 ri— 1 ri— 1 



ni=0 n2=0 n3=0 



+ J{x)6n,i + g d x' Jm{x')D{x' — x)'0„_i(x)r, 



with 



rii + n2 + = n ; rii < n, n2 < n, < n. 



(15) 



(16) 



Notice here that J„ and J„ are given in terms of ilJm{x) and ilJm{x) with m = 0, 1, n — 1 
which are lower order in J . Here 



i?(x, x') = A(x, x') 

- {g'?TU{x) 



j d^X2j d^X3D{x-X3)U{x2)TA-\x2,X3)TU{x3) 



D{x' -X2)U{x')V. (17) 



If ipn{,x) and V^ri(a;) were independent then B{x,x') would be replaced simply with 
A{x,x') and only the first terms survive in Eqs.(0) and ([13|) . 

Once we solve Eq.(|^) for U{x), then we can find i^nix) and ipn{x) for all order n. Now 
the meson field ip{x) can be found using Eq.(|D. Using these results, we can evaluate the 
path integral of Eq.(^ in SPA and find the corresponding action integral S{J'). 

Since only up to second order terms are needed in calculating propagators of each field 
or the energy and the density of a system, let's look at up to the second order in the external 
source J . Dropping the integral sign and the irrelevant ^/-independent terms, the action 
integral becomes up to the relevant second order in J as, 

S{J) = -U{x)J{x) - J{x)U{x) - J{x')B-\x',x)J{x) 
—gU{x)TU{x)D{x — x')Jm{x) + -Jm{x)D{x — x')J„i{x) 
-g'^Jm{x)D{x - xi)U{xi)TB^^{xi, X2)VU{x2)D{x2 - x')Jm{x') 

+ -g'^U {x2)TA~^{x2, xs)TU{xs)D{x3 — Xi)Jm{x4)D{x2 — xijU {xi)T B~^ [xi, x)TU {x)D{x — x)Jm{x') 
-\--g^Jm{x)D{x' — x)U{x)TB^^{x, xi)TU{xi)D{xi — xs) Jm(x4)D(x4 — X2)U{x2)TA^^{x2, X3)TU (xa). (1^ 
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Using Eq.(|^) with W{J') = e*'^^^^ of Eg . ([I8|) , we can find the expectation value of any 
operators up to two body form. 

The mean nucleon field and the mean meson field become 



W{J) \dJ{x) J 
1 ( id 



W{J) \dJ, 




U{x), 



J=0 



(fx' D{x - x')U{x)TU{x'), 



(19) 
(20) 



as expected for a mean field theory. The mean meson field in a nuclear system is generated 
as a virtual meson due to the (real) nucleon field. The meson field propagator becomes, 
dropping the obvious integral sign, 



zA(x — x') = {TLp{x)ip{x')) 



id 



id 



J=o 



W{J) \dJm{x)J \dJm{x')J 

= -iD{x - x') + g'^D{x - xi)U{xi)TU{xi)U{x2)TU{x2)D{x2 - x') 
+ 2ig^D{x - xi)U{xi)TB~\ Xi, X2)TU {X2)D[X2 — x ) 

+ig'^D{x3 — x)U{x2)TA^^{x2, X3)TU{x3)D{x2 — Xi)U{xi)TB~^{xi, X4)TU{x4)D{xi — x') 
+ig^D{x - Xi)U{xi)TB~\xi, xi)TU{xi)D{xi - X3)U{x2)TA-\x2, X3)TU{x3)D{x' - 0:2), (21) 

which has the free meson part and the nucleon field dependent part similarly as in other 
approximations. In RMF, only the second term appears. On the other hand, the propagator 
of nucleon field becomes 



iG{x,x') = {Til){x)il){x 



id 



id 



iB- 



W{J) \dJ{x) 
\x,x') + T[U{x)U{x')]. 



dJ{x') 



W{J) 



(22) 



The second term is the density dependent propagator which is also appearing in RMF or 
RHF. The first term is the correction of the Feynman propagator iA~^{x,x') of the RHF. 
One of the differences of our method from other mean field approaches comes through the 
difference -B(x, x')—A{x^ x') which corresponds to the ladder diagram included in the DBHF. 

The second term of Eq. ([21| ) and the second term of Eq. ([2^ ) form the RHF while the first 
term of Eq.(|22|) gives a correction to the vacuum fluctuation of the RHF. The last three 
terms of Eq.(|2TD gives one and two baryon loops with the Feynman propagator iB~^{x, x') 
for one of the internal baryon line. This propagator contains higher order correlations (such 
as ladder diagram) which were missed in the Feynman propagator 2A~^(x,a;') of the RHF. 
The particle-hole correlations can also be considered in this model through the Feynman 
propagators iA~^{x,x') and iB~^{x,x') in Eqs.(|18|) and (pT|). 

If we neglect the coupling term of '0„(x) and '0„(2;) then there would be no higher 
order correlation effects in the Feynman propagator than RHF, i.e., we would have B{x, x') = 
A{x, x'). If we neglect the nuclear external source, i.e., set J{x) = in Eq.(|^), then the third 
term of Eq.(|r^) does not exist. This means that the vacuum fluctuation comes in through 
the J dependence in a mean field theory. On the other hand, if we neglect the mesonic 
external source, i.e., set Jm{x) = in Eq.(||), then only the first four terms of Eq.(|T8|) 
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survive. This approximation gives the RHF with nuclear vacuum fluctuation. If we neglect 
both J{x) and Jm{x), then we have the RHF without vacuum fluctuation. If we neglect 
any J' dependence of ifj in Eq.(^, i.e., replace ^(a;)r'0(x) in Eq.(^ by lf{x)TU{x), then we 
get the RHF with the flrst flve terms in Eq.(|l8]). However if we use the expectation value 
(u{x)TU{x)'j for ^/'(x)r'0(x) in Eq.@, then we have the RH with vacuum effect. Various 
quantum correlation effects come into a mean fleld approach depending on how the source 
dependence of flelds are treated. 

In conclusion, we have shown that the SPA of the full quantum fleld theory with exter- 
nal sources gives quantum correlation effects which were not included in other mean fleld 
approaches. We also have seen that various assumptions on the external source dependence 
reduce our SPA to RMF, RH, and RHF. Since these further assumptions on the source 
dependence are not perturbative, the various mean fleld approaches of RMF, RH, RHF, 
DBHF, and SPA can not have perturbative relationships. The correlation of ladder dia- 
gram is included in the SPA through the self energy of the propagator iB~^{x, x') without 
any further self-consistency condition in contrast to the DBHF where the ladder diagram is 
included through the self-consistent Bethe-Salpeter equation. Only the HF wave function 
U{x) is required the self-consistent calculation in the SPA in contrast to the DBHF which 
uses the Dyson's equation for the nucleon propagator and the Bethe-Salpeter equation for 
the effective interaction. We should investigate further for the differences between the SPA 
and the DBHF by analyzing more details and by applying the model to a nuclear system 



2T| . To study n-n interaction with this model, we need to consider the action up to 4-th 
order in J'. We also need to extend the SPA further to the case with nonlinear meson fleld 
to consider a nonlinear sigma model or a quark-gluon system. This extension would enable 
the SPA to consider many body interaction within the mean fleld theory. 
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